Selectable linear or quadratic coupling in an optomechanical system 
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There has been much interest recently in the analysis of optomechanical systems incorporating 
dielectric nano- or microspheres inside a cavity field. We analyse here the situation when one of 
the mirrors of the cavity itself is also allowed to move. We reveal that the interplay between the 
two oscillators yields a cross-coupling that results in, e.g., appreciable cooling and squeezing of the 
motion of the sphere, despite its nominal quadratic coupling. We also discuss a simple modification 
that would allow this cross-coupling to be removed at will, thereby yielding a purely quadratic 
coupling for the sphere. 
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Several early optomechanical systems [1-3] were 
straightforward realisations of ideas first proposed by 
Braginsky and Manukin [4, 5], and took the form of a 
Fabry-Perot cavity with one moving mirror. These sys- 
tems are naturally modelled by the Hamiltonian [6, 7] 
(we will use units such that h — 1 throughout this paper) 



H = p 2 /2m- 
thc term 



l x /2+iJii n , where we have introduced 



g\x)a)a 



(1) 



with a the annihilation operator representing the cavity 
field (frequency uj c ), x [p] the position [momentum] oper- 
ator for the moving mirror (frequency ui m and mass m), 
and <7i the optomechanical coupling constant. Physically, 
-ffiin admits a remarkably simple interpretation: the in- 
stantaneous position of the mirror changes the length of 
the cavity, and thus its resonance frequency. 

The light-mirror interaction inherent in Eq. (1) enables 
the passive cooling of the mechanical system, an effect 
that has been demonstrated in a vast array of experimen- 
tal settings [1, 2, 8-10] and that has been progressively 
pushed to the point that optical cooling to the mechanical 
quantum ground state was recently achieved [11]. More- 
over, it has long been recognised that the mechanical and 
optical fields can be prepared in an entangled steady- 
state [12, 13] by the means of Eq. (1). Bhattacharya and 
Meystre [14] proposed a different paradigm, where the 
moving mirror was translucent and placed inside the cav- 
ity, a configuration that was implemented in the seminal 
experiment reported in Refs. [15, 16]. A remarkable fea- 
ture of this model is that, under the right conditions and 
differently from the above, the light-mirror interaction 
may be described through an optomechanical Hamilto- 
nian quadratic in the position operator of the mechanical 
system and reading 
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where gi is a coupling constant. This model has several 
interesting properties, including the possibility of per- 
forming QND measurements of the occupation number 
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of the oscillator [17], but precludes cooling to leading or- 
der in the linearised regime achieved when the cavity field 
is populated by a very large number of photons [18] (see, 
however, Ref. [19] for a treatment that goes beyond the 
linearised dynamics of the system) . 

In this paper we propose to combine these two scenar- 
ios, somewhat in the spirit of proposals that place a cloud 
of ultracold atoms inside an optomechanical setup [20- 
22]. Let us consider the system shown in Fig. 1, which 
includes a small dielectric scatterer (herein referred to 
as the sphere) held in a harmonic potential inside a cav- 
ity [23-29] , one of whose end-mirrors is allowed to oscil- 
late. As we shall show, this system allows us to exert a 
degree of control on the state of the motion of the sphere 
through its indirect interaction with the end-mirror, me- 
diated by the cavity field. As a result of this, and by 
means of a judicious choice of operating parameters, we 
will show that passive cooling of the motion of the sphere, 
despite its apparent quadratic coupling to the cavity field, 
and squeezing of this same motion are both possible. 
Moreover, we expect that the interaction between the 
cavity field, sphere, and end-mirror gives rise to steady 
states with interesting genuinely tripartite entangled [30] 
or nonlocal [31] properties. 

One might think that the Hamiltonian that describes 
this compound system would simply be given by the sum 
of the interaction terms in Hn n through the position x\ of 
the end-mirror and -ff qua d through the position X2 of the 
sphere. As natural as this might seem, this would not 
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FIG. 1. Schematic model of our system: a dielectric nano- or 
microparticle is coupled quadratically to a cavity field that is 
bounded by a moving mirror. The particle is held in place 
by an externally-generated harmonic potential, shown dot- 
ted, whose origin is left open in this work. The cavity field is 
pumped and decays only through the immobile mirror. Dif- 
ferent pumping configurations are explored in Appendix A. 
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be entirely correct, given that the position of the end- 
mirror not only defines the resonance frequency of the 
cavity (as captured through the interaction in H\i n ), but 
also determines the relative position of the cavity field 
(anti)nodes and the sphere. This means that, for the 
configuration shown in Fig. 1, i.e., for a cavity pumped 
from the leftmost (fixed) mirror, the light field couples 
quadratically to the relative position operator x± — x 2 . 
We can thus write the Hamiltonian 



ip 



3i x i 



+ 92 ( 



Xi 



x 2 ) a'a, 



(3) 



where the subscript 'lp' reminds us that the system is 
pumped from the left. In Appendix A, we examine the 
other possibilities for the geometry of the pump field. 

The remainder of this paper is structured as follows. 
Section I will set the stage for exploring the dynamics 
of our system, deriving in particular the Hamiltonian 
and the linearised equations of motion. Section II subse- 
quently makes use of numerical simulations to discuss the 
resulting dynamics. As an application of our scheme, we 
show that the cross-coupling between the two oscillators, 
cf . Eq. (3) , gives rise to efficient cooling of the mechanical 
motion of the sphere. In Sec. Ill we discuss squeezing of 
the motion of the sphere and demonstrate that this mo- 
tion can be squeezed appreciably if the coupling constant 
g 2 is large enough. Finally, we conclude with an outlook 
of possible future directions. 



I. DYNAMICS IN THE LINEARISED REGIME 



where, as before, g\ and g 2 are optomechanical coupling 
constants and \ = xo. 1/^0,2- The signs in front of the 
coupling constants are chosen to aid interpretation: the 
system shown in Fig. 1 has positive g\ , and g 2 > implies 
that the sphere oscillates about an antinode, i.e., a point 
of stable equilibrium. A full description of the system 
requires also the free Hamiltonian, which we write in a 
frame rotating at the frequency of the driving laser wl, 



free 



x 2 



(% 2 
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as well as cavity-pump (i? pU mp) and dissipation Hamil- 
tonians (-ffdiss)j which will be left unspecified. Here, 
A = — w c is the cavity-pump detuning. The full 
Hamiltonian is thus H = iJi nt + -fff rco + H c \i SS + ffpump- 
The corresponding equation of motion for the cavity field 
operator is 

a = (iA-n c )a+i g 1 xi-g 2 (xxi-x 2 ) 2 a-^2n c a in , (8) 

where k c is the cavity amplitude decay rate, and Oi n is 
the input field. As for the sphere and mechanical mode, 
the dynamics is encompassed by 



Pi = -Ui%x + (.91 - xX)a*a - 271P1 - V 2 7i£i- 
P2 = -u 2 x 2 + Xa)d - 272P2 - \f2rf2§2. 



(9) 



Let us start by considering the model in Fig. 1. Clas- 
sically, the intensity of the field within the cavity at 
a distance z from the right mirror is given by I(z) — 
Iq sin 2 (fcz), where the Iq depends on the value of the 
wavenumber k and the amplitude of the input field. To 
treat our oscillators quantum-mechanically we perform 
the replacement 



z — > x\ — X2 + d, 



(4) 



where d is the distance between the equilibrium positions 
of the two oscillators. Let us choose d such that kd — 
nX/2 (n E Z), with A = 2ir/k the wavelength of the 
input light. With this choice, the sphere oscillates about 
a node or an antinode of the cavity field. The interaction 
Hamiltonian involving the sphere thus takes the form 



H cx sin 2 [fc( 



Xi 



x 2 )]"'-- 



a'a. 



(5) 



We now expand Eq. (5) as a series, stopping at the low- 
est order, introduce dimensionless mechanical operators 
[obtained by dividing out the extent of the zero-point 
fluctuations x j = ^h/(rrijU}j) of each oscillator, where 
rrij is the mass and ujj the oscillation frequency of the 
mirror (j — 1) or the sphere (j = 2)], and add back the 
linear optomechanical coupling between a and x\. This 
gives us the total interaction Hamiltonian 



-91X1 
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x 2 ) a'a. 



(6) 



where X = 2g 2 (\xi ~ ^2), 7j are the mechanical decay 
rates, and £j is the self-adjoint zero-mean operator de- 
scribing the Brownian noise affecting oscillator j = 1,2 
and characterised by the correlation functions [32] 



(^(*)6(* , )> = (2«i + i)M(*-*'), 



(10) 



where rij — [ e hu i/( k BTj) _ j s ^ e equilibrium 

phononic population of oscillator j at temperature Tj 
(&B is the Boltzmann constant). 



A. Classical steady-state solution 

Equations (9) are tackled in two steps. First we re- 
place each operator 6 by its quantum-mechanical mean 
value o. The resulting equations of motion are identi- 
cal to the classical ones. Their solution is then used in 
a second step, as described in Sec. IB, to analyse the 
linear evolution of the system. We remark here that 
£j =0 (j = 1,2), and that a m relates to the input 
power as P in = huj L \a in \ 2 . When the dynamics is stable, 
the long-time classical solutions are 



iA 



x 2 = X\- 



2g 2 \a\ 2 2g ig2X \d[ 
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where we have replaced the detuning A with A = A 
giXi — g2\X%i ~ X2) and introduced the quantities 



fi, =ij 1 +2g 2 x 2 \a\ 2 
f2 2 = ^2 + 2g 2 |a| 2 - 



?lx 2 |a| 4 



W2 



2g 2 \ar 



(12) 



Equations (11) diverge in a notable way from the anal- 
ogous ones reported in Ref. [18]. The latter case corre- 
sponds to ours upon going to the infinite-mass limit for 
the end-mirror (\ = gi = 0). From Eqs. (11) it then fol- 
lows that x 2 = Xi = 0. An immediate consequence of this 
is that in the linearised equations of motion, as obtained 
in the next Section, the motion of the mirror is decoupled 
from the rest of the system [cf. Eq. (Bl)]. Under these 
conditions, therefore, the effective mechanical decay rate 
of the sphere cannot be altered by the dynamics in the 
linear approximation, i.e., the motion of the sphere can- 
not be cooled by a quadratic interaction with the cavity 
field. Outside the infinite-mass limit, however, the two 
oscillators interact to avoid this pitfall. As we shall see 
later on, one can choose the parameters such that cool- 
ing of the sphere to a temperature determined mostly by 
the properties of the end-mirror is indeed possible, thus 
demonstrating the control mechanism at the centre of our 
work. 



B. Linearised equations of motion 
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FIG. 2. Occupation number for the motion of the sphere, pre- 
sented on a logarithmic scale. Note the 'resonant' condition, 
indicated by the dashed line, corresponding to the sphere fre- 
quency that minimises this occupation number at each partic- 
ular wi . Cooling of the motion of the sphere by almost three 
orders of magnitude is seen, despite the nominal quadratic 
coupling. The plot here has been optimised over A and Pi n . 
The mirror and sphere are in contact with separate baths at 
temperatures of 50 mK and 1 K, respectively. The rest of the 
parameters used for generating this plot are in Sec. II. In- 
set: Cut of the main figure along the line uji = 10k c (blue, 
lower, curve). The red, upper, curve shows the occupation 
number in the absence of any cooling. The dashed line indi- 
cates the value of U2 for which the lowest occupation number 
is reached. 



Having at hand the explicit form of the operators' 
mean values, we now consider the zero-mean quantum 
fluctuations 56 = — 0. When |o| 2 >• 1 one can truncate 
the equations of motion to first order in the fluctuations, 
yielding the set of linear equations of motion 

^R = AR + R in . (13) 

Here, we have defined the vectors of quadrature fluc- 
tuations R = (Sx,5p,Sxi,5pi,5x 2 ,8p 2 ) , where Sx = 
(8a + 8a i )/V2 and Sp = i(6aJ — 5a) /V2, as well as Ri n = 

-(y/^KcdXin, y / 2K^Sp ia , 0, y/^ii^O, y 7 ^^) with 5x in 

and Sp m that are defined analogously to Sx and 
Sp. The input noise field obeys (5ai n (t)Sai n (t')) = 
(5a in (tySa in (t>)) = and (Sa in (t)5al n (1f)) = 5(t - t'). 
Moreover, we have introduced the drift matrix A, whose 
(lengthy) expression is given in Appendix B. Stability 
of the dynamics encompassed by Eq. (13) is easily stud- 
ied by using the well-known Routh-Hurwitz criterion [33, 
§2.3.2], which (in its most straightforward formulation) 
guarantees the existence of a stable solution if all of the 
eigenvalues of A have negative real part. All of the fol- 
lowing analysis is conducted in this stable regime. 

The linearised dynamics for the operator increments 
can be derived from a quadratic effective Hamiltonian. 
For a fixed set of parameters, the dynamics therefore 
leads any initial Gaussian state to a unique final Gaussian 



steady state. This steady state has zero first moment, 
since we operate in a shifted picture involving only zero- 
mean operator increments, and the matrix of second mo- 
ments V = [(R{t) (g) R(t')) + (R(t) <g> R(t')) T }/2 satisfies 
the Lyapunov equation A V + V A T = D, where 

D5{t-t') = ±[(R in (t)®R in (t')) + (R in (t)®R in (t , )) T ] . 

We remark that the Lyapunov equation admits a com- 
pact solution in in terms of the eigensystem of A [33, 
§5.0.4]. Having at hand the covariance matrix V, we can 
readily analyse the (Gaussian) steady state of the system, 
including bipartite entanglement between any two com- 
ponents of the system, genuine tripartite entanglement, 
or occupation numbers for the two harmonic oscillators. 

II. COOLING OF THE SPHERE 

In this Section, we shall demonstrate the claimed cool- 
ing mechanism for the trapped sphere, and the control 
we can operate on it through the means of the mechan- 
ical mirror. Before specifying the values of the various 
parameters that enter our dynamical equations, let us 
give explicit formulae for g\ and g 2 - The basic theory 
of cavity optomechanics identifies g\ with the first-order 
change in the cavity resonance frequency when the end- 
mirror moves by a distance xq^i [34]. For a cavity of 
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FIG. 3. (a) Normal-mode frequencies, extracted from the 
eigenvalues of A. For low powers, the red (top) curve rep- 
resents the cavity mode, the green (middle, at low powers) 
curve the mirror, and the blue (bottom) curve the sphere, 
(b) Occupation numbers for the sphere (blue; upper curve at 
low powers) and mirror (green; lower). Note the narrow range 
of powers for which the systems hybridise. At this point, the 
two frequencies and occupation numbers are equal. We have 
chosen cji = 10jc c and the optimal values for CJ2 and A. The 
mirror and sphere are both in contact with baths at a temper- 
ature of 1 K. The dashed line denotes the value of P[ n for which 
the occupation number is minimised, (cji = 10k c , 0J2 = 3.4k c , 
A = -27.2kc, gi = 1.0 x 10" 3 K C , Qi = -2.4 x 10~ 10 «: c , 
X = 3.7xl0 -3 ; other parameters as in the body of the paper.) 



length L we can thus write [35] 



,9i 



L V rriiUJi 



(14) 



For g 2 , in the case of spheres that are not too large, we 
use the expression in Ref. [27] 



3V 



£r~l 

e r + 2 



(15) 



where V is the volume of the sphere, V c — jW 2 L is 
the mode volume for a cavity with waist w, k — 2tt/X 



is the wavenumber of the driving field, and A its wave- 
length. The sign of this expression is chosen depending 
on whether the sphere is at a node (— ) or an antinode (+) 
of the cavity field. In the case of a lossless dielectric of 



real refractive index n, 
this expression as 

92 = ±12tt 



Re{e r } 



n 2 + 2 L p(\w) 2 uj 2 



and we can write 



(16) 



where p is the mass density of the material composing 
the sphere. Note, in particular, that g 2 is independent 
of the radius r of the sphere. For silica p = 2650 kg m -3 
and, in the near infrared, n ~ 1.5. We shall use A = 
1064 nm, L — 0.5 cm, k c — 2ir x 50 kHz, r = 0.5 |J.m, 
w = 40 (a.m (which can all be found in Ref. [29] ) , and 
rni = 40 ng [36]. For the mechanical damping constants, 
we take the values 71 = 2ir x 140 Hz [3] and 72 = 2ir x 
0.5mHz. At the two nominal frequencies u)\ = 1 MHz 
and ui 2 = 2tt x 200 kHz, we find gi w 2ir x 36 Hz, 172 ~ 
— 2ir x 10 |xHz at a node of the cavity field, and \ ~ 
2.9 x 10 -2 . In what follows we will not constrain ourselves 
to these particular values of ujj , however, and will scale 
gj and x appropriately. In passing, we note that gi and 
g 2 can be made larger by decreasing the length of the 
cavity or lowering the mechanical oscillation frequencies. 
Moreover, g\ can be increased by lowering the effective 
mass of the moving mirror, while g 2 by using a sphere 
with a larger refractive index or a cavity with a smaller 
waist. Finally, the value of x can be made larger by 
increasing the ratios lo 2 /uii and m 2 jm\. 

We have already mentioned how x — gives rise to 
the steady-state solution x 2 — [18]. In the linearised 
regime, the motion of the sphere is thus decoupled from 
that of the dynamics of the cavity-mirror system. This 
can easily be seen by substituting x = and x 2 = 
into the drift matrix A. A straightforward method to 
explore the implications of the coupling induced by our 
interaction Hamiltonian is thus to look at cooling of the 
motion of the sphere. A sample of such data is shown 



in Fig. 2. The cases (i) uj 2 > uj\ or (ii) uo 2 < 



have 



been excluded from the plot for reasons that will be clar- 
ified later on. We see that, for each u>i, there exists a 
band of frequencies lo 2 for which the steady-state occu- 
pation for the spherical motion is significantly below the 
starting value. Moreover, for a particular a opt = uji/uj 2 
that depends most strongly on the properties of the mir- 
ror (in particular, the constant parameter giy/UJi), the 
motion of the mirror and sphere equilibrate to the same 
occupation number. The data displayed in Fig. 2 show 
the reduction of the occupation number of the sphere by 
almost three orders of magnitude. However, in our nu- 
merical exploration we have seen evidence that this is not 
an upper limit to the efficiency of the process. Indeed, 
decreasing the temperature of the thermal bath to which 
the mirror is coupled will yield a correspondingly lower 
occupation number for the sphere, as a result of dynam- 
ical equilibration of the sphere through the mechanical 
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The basic mechanism through which this cooling pro- 
cess occurs is a hybridisation of the motion of the two 
oscillators. For moderately large values of w 1 /cj 2 , increas- 
ing the input power from zero causes the frequencies of 
the normal modes of the system, given by the imaginary 
parts of the eigenvalues of A, to shift. Thus, for instance, 
the mode that at P m = describes the mirror motion 
shifts to a smaller frequency uj\ [37]. At the point where 
u)% ps ll>2, the two modes hybridise and cooling ensues. 
For lji/u)2 < 1, this hybridisation is no longer possible. 
However, if u>i/uj2 is too large, the system reaches in- 
stability before the criterion a)i ps UJ2 can be satisfied. 
Conversely, if uj\/ui2 is not large enough, the power at 
which d)i ps u)2 is satisfied is too low for efficient cooling. 
Such hybridisation can be seen in Fig. 3, where we plot 
the case ui\ = 10k c , at the optimal values for u>2 and A. 

Another unique feature of our system is that, after 
cooling the sphere, we can switch to coupling light into 
the cavity from the moving mirror. As explained in Ap- 
pendix A, this means that the Hamiltonian that couples 
the motion of the sphere to the optical field is purely 
quadratic. This renders possible, for example, quantum 
non-demolition (QND) measurements of the energy of 
the harmonic oscillator embodied by the sphere. This 
protocol requires that the cavity can be pumped from 
either end-mirror, but is deterministic and otherwise less 
demanding than, e.g., measurement-based schemes for 
achieving the same [38] . 
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III. SQUEEZING OF MECHANICAL MOTION 

The coupling of mechanical elements to a light field, 
both linearly [39] and quadratically [19], can give rise to 
squeezing of the mechanical motion, where the quantum 
nature of the optomechanical interaction acts to suppress 
the noise in one of the quadratures of the mechanical mo- 
tion. To investigate the occurrence of mechanical squeez- 
ing in our setup we set the temperature of the baths cou- 
pled to both oscillators to zero. For the parameters given 
above, there is no observable squeezing in the motion of 
the sphere, and to induce such effects we have found it 
necessary to use a larger value for g 2 , e.g., two orders 
of magnitude larger in the case of the data presented in 
Fig. 4. 

To quantify the amount of squeezing in the steady state 
of the linearised dynamics, we used the figure of merit 



3 2min{(^),(^)}' ^ > 

as plotted in Fig. 4(b) for j — 2. In other words, 
Sj > 1 only when the variance of one of the quadra- 
tures dips below the ^-quantum level that is due to vac- 
uum fluctuations in the ground state. We find a maximal 
S 2 ps 1.2 = 0.7 dB ((pi) ps 0.43) just before the instabil- 
ity threshold, which compares favourably to the amount 
of squeezing obtained in Ref. [19] for pure quadratic cou- 
pling. 



FIG. 4. (a) Quadrature variances, {x 2 } (solid curves) and 
(pj) (dashed curves); the green, upper, curves represent the 
mirror (j — 1), the blue, lower, curves the sphere (j = 2). The 
solid black horizontal line represents the ^-quantum ground- 
state vacuum variance. Note that only the p-quadrature of 
the sphere is squeezed below the vacuum level. The baths for 
both oscillators were held at zero temperature, and <?2 was 
a factor of 100 larger than the parameters in the body of 
the manuscript, (b) Corresponding squeezing for the sphere. 
(oji = 20k c , oj 2 = 10k c , A = -10k c , 0i = 7.2 x 10~ 4 k; c , 
g2 = —8.0 x 10 _9 K C , X ~ 4-5 x 10 -3 ; other parameters as in 
the body of the manuscript.) 

We conclude this Section by noting that the larger 
quadratic coupling strengths necessary for generating 
and observing squeezing can be obtained by making use 
of avoided crossings between pairs of cavity modes [16]. 



IV. CONCLUSIONS 

We have explored a system consisting of an optome- 
chanical cavity with an additional intra-cavity particle 
whose motion is coupled quadratically to the cavity field. 
A direct coupling between the two moving elements arises 
naturally in this system. This coupling can be turned 
off and on at will by choosing which port the cavity 
is driven from, and therefore provides complete control 
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FIG. 5. Different pumping geometries in one dimension, 
(a) For symmetric pumping, a node or antinode always lies 
at the centre of the cavity, (b) For pumping from the immo- 
bile mirror, the nodal structure shifts with the moving mirror, 
(c) For pumping from the moving mirror, the nodal structure 
is fixed with respect to the immobile mirror. 



over whether the interaction with the particle is purely 
quadratic, or whether it has some linear character. As ex- 
amples of the potential of our scheme, we have discussed 
the possibility of cooling the motion of the sphere, de- 
spite its nominal quadratic coupling, as well as squeezing 
of this same motion. Our work opens up interesting per- 
spectives for quantum state-transfer between two oscil- 
lators and the distribution of multipartite entanglement. 
Moreover, the system that we have addressed leaves room 
for asking interesting questions on the nonlocal nature of 
the three-mode state achieved via the coupling mecha- 
nisms addressed here, a goal that can be pursued, for 
instance, through the formal apparatus recently put for- 
ward in Ref. [31]. 
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Appendix A: Alternative pumping geometries 

The function of this Appendix is to elucidate the dif- 
ferences between pumping the cavity from the immo- 
bile mirror, from the moving mirror, or from both ends. 
These three possibilities, which are exhaustive in one- 
dimensional geometries, arc illustrated in Fig. 5. The ge- 
ometry used throughout this paper is shown in Fig. 5(b). 
The most 'natural' case, where the cavity field is sym- 



metric with respect to coordinate inversion, is depicted 
in Fig. 5(a). In both these situations, when the mirror on 
the right moves, the nodal structure shifts with respect 
to the laboratory frame, yielding a term in the Hamil- 
tonian that couples the two oscillators. In the case of 
Fig. 5(c), however, the nodes of the field do not shift 
when the mirror moves, such that there is no longer any 
cross-coupling between the two oscillators. The case of a 
resonant cavity pumped from one end is special: since all 
the input light is transmitted and none is reflected off the 
cavity, both mirrors effectively interact with three run- 
ning waves of nonzero amplitude, two moving away from 
the pump, one towards it, such that the right mirror is a 
time-reversed copy of the left mirror. The entire system 
is therefore invariant under inversion of coordinates and 
time reversal. 

In our notation, the part of the interaction Hamilto- 
nian that involves the sphere can be represented as 



H 



i i — X2) 2 a i a (a) 
[Xi — X2) 2 cJ a (b) 
xia) a 



(c) 



with reference to the respective cases in Fig. 5. The main 
text can be adapted for the remaining cases by setting 
X — > x/2 [case (a)] or x — > [case (c)]. 

This interesting result can easily be seen for case (b) 
(and analogously for the other two cases) by calculating 
the field inside the cavity, normalised to the input field, 
at a distance z from the right mirror. That is 

E{z) = t [e lk ^ L -^ + re lk{L+z) + r V fc(3L " z) + . . . ] 

f ikL 

=: C(L) x £{z) , 

where L is the length of the cavity, k the wavenumber of 
the field, C{L) accounts for the spectral profile of the res- 
onance, and £{z) depends only on the distance from the 
right mirror. For simplicity, the two mirrors are assumed 
to be identical, with (real) reflectivity r and transmis- 
sivity t, and any phase shifts upon reflection or trans- 
mission are absorbed in L; these simplifications are not 
crucial to deriving this result. In the good-cavity limit 
r — > —1, £{z) reduces to the expected sinusoidal profile, 
with (anti)nodes lying at a fixed distance from the right 
mirror. 



Appendix B: Drift matrix 

Here we provide the explicit form of the drift matrix 
A for our problem. In order to aid comparison between 
datasets with different g\ and §2, we introduce the ef- 



fective detuning A = A + g\X\ — (?2(x x i 
resulting drift matrix reads 



X 2 



The 



7 



-k c -A 
A -k c 


[fli - 2.92x(x^i - x 2 )]x [gi - 2g 2 x(xxi - x 2 )]p 

2ff2(x^i - x 2 )x 2g 2 {xx\ - x 2 )p 



-[gi - 2g2x(xxi - x 2 )]p 
[.9i - 2ff 2 x(x^i ~x 2 )]x 


-wi - .92X 2 (x 2 + p 2 ) 


g 2 x{x 2 +P 2 ) 






-252 (x^i - £2)p 


" 





2ff2(x^i - S 2 )^ 





Wi 








-2 7 i 


g 2 x(x 2 +P 2 ) 











UJ 2 





-w 2 - g 2 (x 2 + p 2 ) 


-2 7 2. 



(Bl) 



The phase of the input field can be chosen, as is com- 
monly done, to set p — 0. One can recover the usual 



linear or quadratic optomechanics drift matrices by set- 
ting g 2 — or x = .9i = 0, respectively. 
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